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CHAPTER I 


INTRODUCTION 


A moment method (MM) [2] solution is applied to two different antennas in 
the presence of a grounded dielectric substrate. A grounded dielectric substrate is 
a dielectric substrate of thickness T on top of an infinite ground plane. The first 
antenna is a microstrip antenna consisting of a perfectly conducting rectangular 
patch lying horizontally on top of the dielectric substrate. Figure 1 in Chapter II 
shows the geometry of this antenna. The problems of interest for the microstrip 
antenna are plane wave scattering and the input impedance when the antenna is 
excited by a filament of uniform current. The antenna can be loaded by a lumped 
impedance in a filament of uniform current. The other antenna is a perfectly con- 
ducting vertical wire, or monopole antenna, extending from the ground plane, into 
or through the dielectric substrate. Figure 12 in Chapter III shows the geometry 
of this antenna. The problems of interest are the self impedance of a monopole 
and the mutual impedance between two monopoles. 

The MM solution for the perfectly conducting rectangular patch is begun by 
considering the patch to be excited by an arbitrary source with currents (J*,M‘) 
located in either the dielectric substrate or the free space region above it. The 
surface equivalence theorem is used to replace the patch with an unknown surface 
current, J^. The boundary condition that the total tangential electric field on the 
patch surface must be zero is applied to yield the integral equation for the unknown 
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surface current J'^. The unknown surface current J'^ is expanded into a sum of 
N (known) entire domain expansion functions. The MM solution determines the 
N coefficients in this expansion. For plane wave scattering, the patch is excited 
by a uniform plane wave with an incident angle of The scattered field 

is the sum of the field of each basis function multiplied by the strength of that 
basis function. To find the input impedance, the patch is excited by a vertical 
filament of uniform current extending from the ground plane to the patch. The 
input impedance is determined from the ratio of the voltage at the feed point and 
the current into the feed point. The effect of the load is determined from modeling 
the microstrip antenna as a two port device, where the excitation represents port 
1 and the load represents port 2. The current through the load is determined from 
the unloaded case. The load current acts as a source and is added to the unloaded 
excitation to obtain the loaded result. 

Many papers have been written on the rectangular patch microstrip antenna. 
Early work by Munson [3] and Derneryd [4] modeled the microstrip antenna as a 
transmission line that radiates from the open ends. This model does not include 
higher order modes and thus is valid only near the first resonance of the patch. 
Lo et ad. [6] and Richards et al. [7] used a cavity model. This model assumes 
an admittance around the open edges of the patch to form a cavity. Higher or- 
der modes are supported by a cavity model. Agrawal and Bailey [5] performed 
a wire grid analysis of the microstrip patch. This was one of first solutions that 
attempted to model the microstrip patch by the currents on the patch. Newman 
and Tulyathan [8] employed a MM solution similar to that for plates. Bailey and 
Deshpande [9,10] and Pozar [15] used the Greens function for the grounded dielec- 
tric substrate in their MM solutions and used entire domain mode functions. Pozar 
compared the numerical efficiency of the entire domain modes and the piecewise 
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sinusoidal modes when solving for the impedance matrix. Most of the work up 
to this point was primarily concerned with the behavior of the microstrip patch 
in a narrow band around first resonance. Newman and Forrai [14] examined the 
wide band scattering of the microstrip patch. An impedance matrix interpolation 
scheme was employed to minimize the number of times the MM impedance matrix 
had to be computed. This report will examine the wide band scattering and the 
wide band input impedance of a loaded rectangular patch microstrip antenna. It 
will be shown that when the antenna is loaded by a vertical filament of uniform 
current that contains a lumped impedance, the results closely approximate the 
behavior of a microstrip antenna loaded through a coaxial port. 

The MM solution for the perfectly conducting vertical thin wire is begun by 
considering the wire to be excited by an arbitrary source with currents (J*,M*). 
The wire is replaced by a unknown surface current J^. Thin wire approximations 
are applied to this surface current. It is modeled as a vertical or z directed filament 
of current [17]. This current is approximated by a sum of N (known) basis 
functions. The MM solution determines the unknown coefficients in this expansion 
for excitation by the so called delta gap generator [19]. The input impedance is 
determined from the ratio of the excitation voltage and the current on the wire at 
the point of excitation. The mutual impedance is determined from the ratio of the 
excitation voltage on one wire to the current on the other wire. 

Collin and Zucker [18] analyzed a vertical electric current element inside a two 
layer earth. They modeled the earth as flat and therefore their analysis is similar to 
that of a vertical current element inside of a grounded dielectric substrate. Pinhas 
and Shtrikman [22] analyzed vertical currents in the dielectric region of a grounded 
dielectric substrate by modeling the vertical current as a horizontal surface current. 
One of their primary Interests with this analysis was modeling the vertical current 


3 



for a coaxially fed luicrostrip antenna. Chi and Alexopoulos [21] performed an 
analysis of a vertical wire in a grounded dielectric substrate which is similar to the 
one presented in this report. Their excitation was a magnetic frill current in the 
aperture of the coaxial feed on the ground plane. This report will use the simpler 
so called delta gap generator for the excitation and compare the results to theirs. 
It is shown that the two different excitations produce nearly identical numerical 
results. 

Chapter II describes the derivation of the integral equation and the MM solu- 
tion for the microstrip antenna. Chapter III describes the same derivations for the 
vertical monopole antenna in the grounded dielectric substrate. Chapter IV con- 
tains the results for plane wave scattering and input impedance of an unloaded and 
loaded microstrip antenna. Data from both a MM computation and measurement 
is presented. Chapter IV also shows the mutual impedance between two vertical 
monopoles in a grounded dielectric substrate. Chapter V summarizes this report. 
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CHAPTER II 


THE INTEGRAL EQUATION AND MOMENT METHOD 
SOLUTION FOR THE RECTANGULAR PATCH MICROSTRIP 

ANTENNA 


2.1 Introduction 

This chapter will describe the integral equation and the moment method [2] 
(MM) solution to the problem of plane wave scattering and input impedance of 
a rectangular patch microstrip antenna on an infinite microstrip slab. The basic 
geometry of this structure is shown in Figure 1. The antenna can be excited by 
either a filament of uniform current extending vertically through the dielectric 
substrate or by an incident plane wave. A load will be modeled by an uniform 
vertical current with fixed impedance. 

The solution is begun by using the equivalence theorem to replace the perfectly 
conducting patch by an unknown surface current J"^. Then, the integral equation 
for is obtained by requiring that the total electric field tangent to the patch 
surface be equal to zero. Finally, the integral equation is solved using the MM. 

2.2 Derivation of the Integral Equation 

This section describes the use of the surface equivalence theorem to derive 
the integral equation for the problem of radiation or scattering by a perfectly 
conducting rectangular patch on an infinite dielectric substrate on top of an infinite 
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a) Top view of a perfectly conducting rectangular patch on top of a infinite dielec- 
tric slab. 







b) Side view of a perfectly conducting rectangular patch on an infinite dielectric 
substrate on a perfectly conducting ground plane. It is excited by an arbitrary 
source with electric and magnetic currents (J*,M^). 

Figure 1: The basic geometry of a microstrip patch antenna. 
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perfectly conducting ground plane. The geometry of a dielectric substrate on top of 
a perfect conductor will be referred to as a grounded dielectric substrate. Figure 1 
shows a perfectly conducting rectangular patch of length 2H in the x direction 
and width 2W in the y direction on a grounded dielectric substrate. The region 
above the dielectric substrate, referred to as Region J, is free space. The region 
inside the dielectric substrate, referred to as Region II, has a permittivity of C 2 , 
permeability of fiQ, a,nd thickness T. The origin of the coordinate system is placed 
at the center of the rectangular patch. The normal to the rectangular patch surface 
is the z direction. An arbitrary source, which can be located in either Region I or 
Region II, has electric and magnetic currents (J*,M*) that illuminate the patch 
with incident fields (E*,H*). Note that (E*, H*) are the fields of (J*,M*) in the 
presence of the grounded dielectric substrate without the rectangular patch. Let 
(E,H) denote the electric and magnetic fields produced by (J^,M*) in the presence 
of the grounded dielectric substrate with the rectangular patch. 

Since the patch is a perfect conductor, the equivalent current on the rectan- 
gular patch surface will be 

= z X H. (2.1) 

Since the patch is infinitesimally thin, the surface current actually represents 
the vector sum of the surface current on the top surface and of the surface current 
on the bottom surface of the patch [12]. The equivalent problem is illustrated in 
Figure 2. The total electric and magnetic fields in either region are 


E = E* -f- E^ 

(2.2) 

H = H* -1- 

(2.3) 


where (E"^,!!^) are the fields radiated by the equivalent current in the presence 
of the grounded dielectric substrate. 
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a) Top view of the equivalent surface current replacing the perfectly conducting 
rectangular patch. 



: ;K fei gr dJ JiJ W ::: CC a: ; 

b) Side view of the equivalent surface current J'^ replacing the perfectly conducting 
rectangular patch on a grounded dielectric substrate. It is in the presence source 
currents (J*,M*). 

Figure 2: The equivalent geometry of a microstrip patch antenna. 
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On the surface of the patch, the tangential components of the total electric 
held must equal zero. Then, from Equation (2.2) 

(E* + E^) X z = 0. (2.4) 

This can also be written as 

- = EU. (2.5) 

This is the integral equation for since 

E* = /_" dz' dy' (2.6) 

where = G^(x',y’ \ ®,y,z) is the dyadic Green’s function for the grounded 
dielectric substrate. E^ is a function of (®,y, 2 ), while is a function of (®',y'). 
The dyadic Green’s function 6^ is used only to show the integral form of an electric 
field and will not be used in the analysis. 

2.3 The General Moment Method Solution 

This section will use the integral equation developed in the preceding section to 
develop a general MM solution for the geometry of Figure 1. The unknown surface 
current in the integral equation will be replaced by a sum of N known vector 
expansion functions of unknown strengths. By taking an inner product between 
Equation (2.5) and a set of N vector weighting functions, a matrix equation is 
formed. A Galerkin solution will be used and expressions for the matrix elements 
will be given. 

To begin the solution for the unknown equivalent surface current, is rep- 
resented as a sum of N known vector expansion functions 

= E (2.7) 

11=1 
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where are the expansion functions and In represents the unknown strength of 
each expansion function. Let be the field radiated by the expansion function 
in the presence of the grounded dielectric substrate. Then Equation (2.5) can 
be written as 

- ^ ^n^ntan = E|on- (2-8) 

By taking the inner product of Equation (2.8) with a set of N known vector 
weighting functions defined on the surface of the patch, Equation (2.8) reduces to 
a AT X N system of simultaneous linear equations. These equations can be written 
compactly in matrix form as 

IZ][I] = [V] (2.9) 


where [Z] is an JV x JV impedance matrix, [/] is a JV element current column array 
containing the unknown coefficients In of the expansion of J*^ and [F] is a N 
element voltage vector column array. Applying a Galerkin solution, which defines 
the vector weighting functions to be the same as the vector expansion functions, 
a typical element of [Z] is 


fW fH ^ ^ 

Zmn — — E„ ■ 3^ dx dy 

and a typical element of [V] is 

fW [H . . 


( 2 . 10 ) 


( 2 . 11 ) 


By applying reciprocity to Equation (2.11), the elements of the voltage vector can 
be expressed in terms of the electric field of in the presence of the grounded 
dielectric substrate. The reciprocal expression is 


Vm = 


( 2 . 12 ) 
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where the integral is over the volume of the impressed current and is the 
electric field of J^. Now both the impedance matrix elements and the voltage 
vector elements are expressed in terms of the electric field of the expansion and 
weighting functions. 

Note that in the preceding equations, the tan subscripts that appeared in 
Equation (2.8), from which the above were based, have been dropped. This can 
be done since the vector weighting functions have direction tangential to the patch 
surface. 

2.4 Evaluation of the Impedance Matrix 


In this section, exact integral expressions will be developed for the imped- 
ance matrix elements using the general expression developed in the preceding sec- 
tion. The expansion functions representing the unknown surface current will be 
restricted to having either x or y vector direction. Figure 3 shows an expansion 
function, J^, on the interface of Region I and Region II. The electric field for 
J„ in the presence of the grounded dielectric substrate will be developed using the 
plane wave expansion method [13]. Note that this electric field is used for evaluat- 
ing both Equation (2.10) and Equation (2.12). Equation (2.10) will be evaluated 
to obtain the expressions for the impedance matrix elements. 

Since the surface current exists on the interface of Region I and Region II, 
both regions can be considered source free. The fields in each region can be con- 
structed from two scalar functions [Ij. These scalar functions must satisfy the 
scalar wave equation. For Region /, 


(v^ + 4 ) 




= 0 


(2.13) 
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^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 


Figure 3: The MM expansion function, J^, at the dielectric interface. 


and for Region II, 


where, 


+ 4 ) 


^mll 

^tll 


= 0 


«0 = 

«2 = ^^62^*0- 


(2.14) 


(2.15) 

(2.16) 


generates fields with a magnetic field tangent to the planar interface (TM) 
and ^'e generates fields with an electric field tangent to the planar interface (TE). 
The E and H fields of either region can be defined from the ^e, e, and fi of 
the appropriate region using 


E = -V X (i^c) + —V X V X (z^rn) 

JU)€ 

H = V X (z^m) + —V X V X (z’i'c). 

JUffl 


(2.17) 

(2.18) 
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Using the plane wave expansion method [13], the scalar wave functions can 
be defined by their spectral functions. The equations to transform the spectral 
functions to the scalar wave functions are, in Region /, 



1 f°° 

^ml 

^el 

8 

8 

1 

|cs 

II 

1 

1 


dK^dUy 


(2.19) 


where 

- k\- Re(/c^i)> 0, Im(K 2 i)< 0 (2.20) 


and in Region II, 


^m// 

1 

^mll cos(«z2[2 + ^’1) 

^ell 

47T^ J—oo J—oo 

$e//sin(/c^2[z + 


( 2 . 21 ) 


where 

Kz2 — \/«2 - - «y Re(/C22)> 0, Im(/C22)< 0. (2.22) 


The conditions already satisfied are: the boundary condition on the ground 
plane {z x E// = 0) in Equation (2.21), the radiation condition as r — > oo in 
both Equation (2.19) and Equation (2.21), and the integrals in Equation (2.19) 
converge as z — > cx). The following boundary conditions at the dielectric interface, 
z = 0, still need to be satisfied 


z X (E// — E/) = 0 

(2.23) 

i X (H„ - H,) = 

(2.24) 


is assumed to be Fourier transformable and expressed in x and y compo- 
nents. The Fourier transform pair is 


J_ — xJnx "I" y«^i 


ny 
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= 4 ^ '''‘9 ( 2 - 25 ) 


jf = xj+, + yj+y 

= /_^ l”jf [^'^nx + y/ny] dx dy 


(2.26) 


where and J^y are the transforms of Jnx and Jny respectively. 

By substituting Equations (2.19) - (2.22) into Equations (2.17) - (2.18), the 
integral expressions for Ej, H/, E/j, and H// are obtained. Applying the boundary 
conditions of Equations (2.23) - (2.24), with Equation (2.25) substituted into 
Equation (2.24), solves for ^ell terms of Jnx and Jny 

From this result, the total electric and magnetic fields of Region I and Region II 
[13] can be written in terms of Jnx and Jny using Equations (2.17) - (2.18). 

The tangentied electric field components at the patch surface for a x directed 
surface current is 


Exn 

Eyn 


fir 


=-— r 

47T^ J—OO J — OO 


KxKy (Ag — i4j|f) 


and for a y directed surface current 


E, 


xn 


^yn 


4ir^ 7-00 J— 


OO fOO 

oo 


KxKy {Ae — A^) 

{k^ae + r^A^) 


Ji+e-Ji'^^^+'^yyUKxdKy, 

(2.27) 

JS^e-Ji'^-’^-^'^yyUKxdKy 


where 




(2.28) 

Ae 

wfiQ sin Kz 2 T 
k^De 

(2.29) 


KziKz2sinKz2T 

IjJK^Dm 

(2.30) 


_ 2,2 
— Kj. 4* /Cy 

(2.31) 

De 

= Kz2 COS Kz2T + jKzi sin Kz2T 

(2.32) 

Em 

= e 2 K 2 i cos Kz 2 T + }iQiiz2 sin Kz 2 T. 

(2.33) 
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The elements in [Z] are obtained by inserting Equations (2.27) and (2.28) 
into Equation (2.10) [13]. If mode n and mode m are x polarized, the mutual 
impedance is 

^mn ~ I oo I oo ^mx^nx (2.34) 

If mode n is X polarized and mode m is y polarized, the mutual impedance is 

j /OO /OO « , 

^mn — J ooJ ^ ~ (2.35) 

If mode n is 5^ polarized and mode m is 5c polarized, the mutual impedance is 

Zmn — J J i^xf^yl-^E — A^] dKx d.Ky (2.36) 

If mode n and mode m are f polarized, the mutual impedance is 

Zmn = I oo I oo ^ny^tny ^^y (2.37) 

where, for example, is the Fourier transform of J«* (see Equations (2.25) 

and (2.26)) and 

j.-. = dx dy. (2.38) 

Equation (2.39) shows the general form of an element of [Z] and a trans- 
formation of the integral from (Kx,Ky) coordinates to (K,a) coordinates. The 
transformation shown converts two integrals with infinite bounds to one integral 
with finite bounds and another integral with one infinite bound. This simplifies the 
application of the [Z] element expressions since they are integrated numerically. 

/ oo fOO . , 

-oo J—oo 

f2v foo .. . 

= F(/c,a)J^J„ Kdnda (2.39) 

where, 

Kx = K cos O' 

Ky = /csina. 
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(2.40) 

(2.41) 



2.5 Evaluation of the Voltage Vector 


This section will obtain the general equations of the voltage vector of Equa- 
tion (2.11) for an incident uniform plane wave in Region I and for a vertical 
filament of uniform current in Region II. 

2.5.1 Voltage Vector for an Incident Plane Wave in Region I 

This subsection will evaluate Equation (2.11) for the case of a uniform plane 
wave. Figure 4 shows the side view of a plane wave incident on a microstrip patch. 
A nonphysical impressed current that radiates a uniform plane wave in free space 
will be defined so that the voltage vector elements can be written in terms of 
the far zone electric field of a weighting mode in the presence of the grounded 
dielectric substrate. Finally, the far zone electric field of the weighting mode will 
be evaluated asymptotically so that the voltage vector elements can be computed 
analytically. 

The electric field of a uniform incident plane wave in free space can be written 
as 

E^ = Eoe^*' (2.42) 

where 

K = XKQ sin ^ cos -h y«o sin ^ sin -f z/CQ cos ^ (2.43) 

r = ilr sin ^ cos ^ H- yr sin ^ sin ^ -F zr cos ^ (2.44) 

with (f> defining the spherical angle of incidence. An infinitesimal dipole source in 
free space with current of amplitude will produce the plane wave described 

by Equation (2.42). Note that this source is physically unrealizable since the 
amplitude of the source is proportional to the distance r. Since an impressed 
current for a uniform plane wave in free space has been defined and the electric 


16 


Re«g'i_orx I 




^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 


i—rc 


Figure 4: Plane wave incidence on a microstrip patch. 

field of in the presence of the grounded dielectric substrate was developed in 
the evaluation of the impedance matrix elements, Equation (2.12) can be used to 
express the elements of the voltage vector for a uniform plane wave in free space 
incident on the microstrip patch. Defining the impressed current to be in the far 
zone of the microstrip patch, let the far zone field of be 6^. Then, evaluating 
Equation (2.12), the voltage vector elements are written as 


47T 

Vm = e^-Eo. 

JWflQ 


(2.45) 


Using the method of stationary phase [13], the far zone electric field of mode 
m can be evaluated asymptotically. With this evaluation, Cm can be written as 

e-JKo»’ 

— (2.46) 


Bm = [cos (ffCy — sin ^ej.] cos 0 + 0 [cos <f>ey + sin | 

where, for x polarized modes. 


ex = -C[{€ri4-KlKzi cos Kz2T + 3{kI-kI)k^ 2 sin k^2T]J^^x (2-47) 


(2.48) 


iy — CK3.Ky\iiz\cosKz2T jsmK;^2T]Jmx 

and for y polarized modes, 

e® = CKxKy [k^i cos K 22 T + J sin ^my (2.49) 

iy = -C* [(cr«o - ‘^osac 227’ -f j(ko - «y)«z2 sin /C22^| (2.50) 

with 

^ ^ j sin Kz 2 T 

When evaluating Equations (2.47) - (2.50), it is necessary to replace («x,Ky) by 
their values at the stationary phase point [14], 

(«®, «y) — > (/Co sin 0 COS <f), KQ sin 0 sin (f)) (2.52) 

2.5.2 Voltage Vector for a Vertical Uniform Filament of Current in 
Region II 

This subsection will obtain the voltage vector elements for a uniform vertical 
current filament. Figure 5 shows a vertical filament of uniform current exciting the 
patch. Since the impressed current of this source is known and the electric field 
of in the presence of the grounded dielectric substrate has been developed, the 
voltage vector elements will be obtained using Equation (2.12). For this problem, 
the electric field of can not be simplified by a far zone asymptotic approxima- 
tion. The voltage vector element for this excitation is obtained in integral form. 



Starting with Equation (2.12) and letting 


J* 


zIq6{x - X f)S{y - y f) -T<z<0 
0 elsewhere 


(2.53) 
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a) Top view of a microstrip patch excited by a current filament. 
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b) Side view of a microstrip patch excited by a current filament. 

Figure 5: A rectangular patch microstrip antenna excited by a vertical current 

filament. 
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the voltage vector element becomes 


/■O 

^rn = J ^Ezml^dz. 


(2.54) 


Ezm is 


Ezm=-~^ f f COS Kz 2 {z + T)dKxdKy (2.55) 

J-ooJ-cx> uDm “ 


Ky cos + ^) ^Ka, dKy. (2.56) 


for X polarized modes and 
P 1_ f°° f°° 

47 T^ y ~ 00 /-OO UfDj^ 

for y polarized modes [13]. 

Substituting Equations (2.55) - (2.56) into Equation (2.54) and integrating 
over z defines the integral equation for the voltage vector element of a vertical 
uniform current. For x polarized modes, 


1 foo roo Kzi sin K 22 T 


1 roo yc 

47 T^ J—oo J— 


-00 J—OO 

and for y polarized modes, 

foo roo sin Kz 2 T 


dKz, dK, 


1 f°o rc 

= “J? /-«. L 


n * dK^ dKy 

00 l<}K22Dm 


(2.57) 


(2.58) 


Equations (2.57) - (2.58) contain an exponential term dependant on the (x,y) 
location of the vertical current. This factor will produce an oscillatory conver- 
gence in the integration which will complicate its computation with respect to the 
computation of the impedance matrix elements. 

Equation (2.59) shows the general form of the voltage vector element for the 
vertical uniform current. It also shows the transformation from {Kx,Hy) to (k., a) 
coordinates used in the evaluation of the integral. 


J—OO J—oo 


-00 ^—00 
r2r too 


^ lo lo KdKda 


(2.59) 
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The Vfn for a wire of radius a is easily developed from the solution for Vm of 
the filament if the wire is modeled as an infinite number of filaments as radius a 
around some point (a5c>yc)* Using Equation (2.59), the voltage vector of the wire 
is 

2ir Ja Jo Jo 

(2.60) 

By evaluating the ^ integral, Equation (2.60) becomes 

Vrn = f f Jo(«o)«dKda (2.61) 

where Jq is the zero order Bessel function. 

2.6 Loading of a Rectangular Microstrip Patch 

This section will describe the procedure for modifying the voltage vector when 
a lumped load at an arbitrary position («£, y^) is attached between the rectangular 
microstrip patch and the ground plane. The load is assumed to be in a vertical 
filament from z = —T to z = 0 with a uniform current of magnitude Ii flowing 
through it from the patch to the ground plane. The load current is determined by 
modeling the microstrip antenna as a two port device. 

Equation (2.9) describes the matrix equation for the MM solution of the un- 
loaded microstrip patch. When a load Zi is added. Equation (2.9) becomes 

[Z] [/) = [r] - h [v‘‘] ( 2 . 62 ) 

where [v*] is the excitation voltage vector, Ii is the magnitude of the uniform 
current through the load, and is the voltage vector for a 1 amp current 

filament at the load position, (®£,yx;). Note that the 1 amp current from which 
is obtained is in direction opposite to the load current. 
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Figure 6: The equivalent circuit of a symmetric two port with an I{ current 
source at port 1 and a load at port 2. 

In order to determine Ii, the microstrip antenna is modeled as a two port 
device. Figure 6 shows the equivalent circuit of a symmetric {Z12 = Z21) two port 
device excited at port 1 by an Ii current source and with a load of impedance Zi 
at port 2. From circuit theory, the current through Zi is 

" Z22 + 

Since the product of Z12 and /j is the voltage at port 2 with an open circuit at 

port 2, Equation (2.63) can be written as 

I Vn 

^ Z22 + Zl 

where V12 is the open circuit voltage at port 2. 

For the microstrip patch, V12 is equivalent to the voltage at {xi,yi) for an 
unloaded antenna. This voltage is called and can be found from 

Vl = - E, dz (2.65) 
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where Ez is the z directed electric field at This integral can also be 

expressed as 

Vf = - 5; /„ /" dz. (2,66) 

m=l 

where Ezm is the z directed field of patch mode m at (xi^yi) and Im is the 
strength of mode m for the unloaded antenna. Since the elements of are 

obtained for a 1 amp z directed current filament at {xi,yi), Equation (2.12) can 
be written for [V’^j as 

• z dz 

= J ^ Ezm (2.67) 

where Em is the total field of mode m at {xi, yi). From Equations (2.66) and (2.67), 
the matrix equation 

Vr = -W* \y^] (2-68) 

can be written where [/] is the current column solution of Equation (2.9) for the 
unloaded antenna. 

Z 22 for the microstrip patch is obtained by solving for [/] when the excitation 
voltage vector in Equation (2.9) is [F^j. Using this [/] in Equation (2.68) to solve 
for the voltage at the load port, Z22 is the ratio of the voltage at the load port 
to the input current at the load port (which is 1 amp). Note that the excitation 
voltage vector is not used in determining ^22- 

2.7 Expansion and Weighting Functions Used in the MM Solution 

The expansion and weighting functions used in this MM analysis will now be 
defined. The weighting functions are chosen the same as the expansion functions, 
and thus a Galerkin solution is employed. The expansion functions chosen will 
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be entire domain, which implies that each mode exists over the entire rectangular 
patch surface. Each expansion function will be separable into a function of x and 
a function of y. These functions of x and y are Fourier transformable, as required 
in the derivation of the field equations for a surface current on the patch. Below 
the expansion functions and their transforms are presented. The advantages of 
choosing these expansion functions will also be discussed. 

Starting with the x polarized expansion function, a mode n is defined to be 


= X^(x,H)Y,(y,W) 

= An ^(x + H) cos ^(y+W) (2.69) 

where now mode n actually specifies an integer doublet pq where p = 1,2,3,... 
and g = 0, 1,2, ... are the bounds. When the Fourier transform, defined by Equa- 
tion (2.26), is applied to Equation (2.69) the result is 


where 




Y± ^ COS KxH 

^ ~ 4- 

sin KxH 

-J. _ , j2Ky cos KyW 

’ ■ -i-mf 

2/Cy sinKyW 


Kt 


(^) 


Similarily, the y polarized modes are 


P= 1,3,5,. 
p = 2,4,6,. 
q = 1,3,5,, 
g = 0,2,4,. 


•frny — X^{y,W)Y^,{x,H) 




(2.70) 


(2.71) 


(2.72) 


(2.73) 
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with the corresponding Fourier transform 


W)Yf(K.,,B). (2.74) 

Note that p' ,q' can be the same or different from p,q. 

The transforms of the expansion and weighting functions are composed of 
even and odd functions of and Ky. Realizing that the integral from — oo to oo 
of an odd function is zero, certain impedance matrix and voltage vector elements 
may be identically zero. As an example, consider the impedance matrix integral 
for an element where both modes n and m are x polarized. For this integral, the 
function F{Kx,Hy) (refer to Equations (2.34) and (2.39)) is an even function of k® 
and Ky. Using Equations (2.71) - (2.72), Figure 7 shows whether the expansion 
function Jnx >s even or odd with respect to Kx and Ky for all possible combinations 
of p and q being even or odd. Figure 8 shows whether the integrand is even or 
odd with respect to k* and Ky for all possible even-odd combinations o{pn,qn,pm, 
and qm and whether the integral is non- zero or not. Only 1/4 of the elements in 
Figure 8 are non-zero. This allows a simplification of the general integral described 
by Equation (2.39). For the non-zero elements, the a integration is performed from 
0 to 7t/2 and the result multiplied by 4. The zero elements are set to zero. Similar 
results are found with Equations (2.35)- (2.37). 

With the voltage vector integrals of Equations (2.57) - (2.58), the exponential 
term can be separated into functions of k® and Ky and these functions can be 
broken into a sum of an even function and an odd function. Equation (2.75) shows 
this process. 

^-]{KxXf+Kyyf) _ ^-JKxXf^-JKyVf 

= {cosKxXf — jsinKxXf){cosKyyf — jsinKyyf). (2.75) 
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Figure 7: The /c® and Ky functions of Jnx as even or odd functions with respect 

to p and q. Note: e is even, o is odd. 

Expanding Equation (2.75), the integrand can be written as the sum of four terms, 
each separable into functions of and Ky. By analyzing the even-odd functions 
in each of these terms, only one term is non-zero. Figure 9 shows the integrands of 
an X polarized mode voltage vector for all possible p and q even-odd combinations. 
For the voltage vector, none of the elements are zero in general, however, by having 
integrands that are even functions of and Ky, the a integration can be performed 
from 0 to ir/2 and multiplied by four. 

One advantage of the entire domain expansion functions is that as the patch 
increases in electrical size, the width of the current transforms decreases. This 
results in more rapid convergence of the integrals involving the current trans- 
forms. The mode shapes also correspond to the resonant current distributions on 
the patch. This enables a simple physical interpretation of the magnitude of the 
elements in the current column vector. 

2.8 Computation of the Scattered Field and Input Impedance 

Once the impedance matrix elements and voltage vector elements are known, 
the current column array can be obtained using standard matrix algebra. Assum- 
ing the current column is known, this section will obtain expressions for the far 
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Figure 8: The and Ky functions of the integrand (F(«x>Ky) is even for both 
Kx and Ky) as even or odd function with respect to pn^gniPm and qm- 
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Figure 9: The integrands for a x polarized mode voltage vector with respect to p ■ 

and q. | 

zone fields of the patch. These far zone fields are referred to as the radiated fields | 

when the excitation is by the vertical current filament and as scattered fields when 
the excitation is by an incident plane wave. For excitation by the vertical uniform 
current, the input impedance of a mictrostrip antenna will be found. 

The electric field of the patch is found from the sum of the field of each mode 
weighted by the strength of each mode, 

N 

E = 5 : /„E„ (2.76) 

n=l 

where En is the radiated field of mode n as given by Equations (2.27) and (2.28). 

For the far zone electric field, Equation (2.76) reduces to 

N 

e(<?,<^)= E (2.77) 

n=l 

where Cn is the far zone asymptotically approximated field of mode n, defined in 
Subsection 2.5.1. 

The input impedance of a microstrip antenna is the ratio of the voltage at the 
feed port to the input current. The voltage at a port on the microstrip patch for 
arbitrary excitation of the antenna is given by Equation (2.68) in the Section 2.6. 

When the excitation is a 1 amp uniform current at the port where the voltage is I 
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desired, the input impedance can be written directly from Equation (2.68) as 

= -[I]^[V] (2.78) 

where [/] is the current column solution and [V] is the voltage vector at the input 
port. The probe itself tends to be inductive, and the reactance of this inductance 
can be modeled as [9] 

= (2.79) 

where /mHz frequency in megahertz. Using Equations (2.78) - (2.79), the 

input impedance becomes 

Zin=-[lf[V\+jXi. (2.80) 

The input impedance of a microstrip patch excited by a vertical filament of 

uniform current shows good agreement to the input impedance of a coaxially fed 
microstrip antenna when the thickness of the substrate, T, is much smaller than 
one wavelength. 

2.9 Numerical Difficulties 

There are several numerical problems to overcome when computing the plane 
wave scattering or input impedance of a rectangular microstrip patch antenna. 
One problem is the CPU time to compute the impedance matrix. This problem 
is especially severe if one desires to compute the wideband scattering or input im- 
pedance of the antenna. Both the impedance matrix and voltage vector element 
integrands contain surface wave poles which produce singularities in the integra- 
tion. Finally, there is the exponential term in the voltage vector element integral 
which produces an oscillatory behavior of that integral. Each of these problems 
will be discussed in the following subsections. 
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2.9.1 CPU Time Reduction for the Impedance Matrix Computation 

Equation (2.39) shows the general form of an impedance matrix element in- 
tegral expression and it is seen that the integrand is a separable function of 
and Jn- To generate all the elements of the impedance matrix, for every cal- 
culation there must be N calculations of J^. Since there are N calculations of 
the total time to compute all the elements would be proportional to N^. The 
impedance matrix is symmetric, therefore only N^/2 elements need be calculated 
to generate all the elements. 

Referring to Equations (2.71) - (2.72), the transforms and for the cho- 
sen expansion and weighting functions will either be equal or equal and opposite, 
depending on n. This indicates that when is calculated, is also known. 
When computing the elements of the impedance matrix, all N are calculated 
(which means that all N are obtained as well) within the integral expression of 
Equation (2.39). The impedance matrix elements are then obtained by multiplying 
and together for n = l,...,iV and m = l,...,iV. Since the time required 
to multiply and together (which is proportional to N^/2 when taking the 
symmetry into account) is significantly less than the time to compute andJ^ 
(which is proportional to W), the total time to compute all the elements of the 
impedance matrix tends to be proportional to iV. 

2.9.2 CPU Time for Wide Frequency Band Calculations 

To make a broadband computation of the scattering or input impedance of a 
microstrip antenna, many frequency points are needed to show the behavior near 
resonance since microstrip antennas are narrow band. Thus, although the CPU 
time to evaluate the matrix elements at a single frequency may be tractable, the 


30 


CPU time required to evaluate the matrix elements at hundreds of frequencies over 
a wide bandwidth can be prohibitive. 

A way to make broadband calculations of a microstrip antenna starts by an- 
alyzing the impedance matrix with respect to frequency. The elements in the 
impedance matrix are very slowly varying with frequency [16], which is quite dif- 
ferent from the scattered field or input impedance data. An interpolation tech- 
nique can be employed to approximate the impedance matrix elements at fre- 
quencies in between the frequencies at which the elements are calculated using 
Equations (2.34) - (2.37). Over the desired bandwidth, data points separated by 
a frequency interval of Afz are required to characterize the behavior of the micro- 
strip antenna. The number of computations required to get data points every Afz 
would be impractical, so the impedance matrix elements are computed at intervals 
of Afc instead. Typically, Afc Afz- Three of the points where the elements 
are calculated directly, /<-i, /c2» /c3> (with fd < fc2 < /cs) will define a frequency 
band of width 2 A fc- A quadratic interpolation method can approximate the val- 
ues of the elements every Afz in between fd, fc2 /c2> /c3- By dividing the 
desired frequency band into sub-bands of width 2Afc and applying the quadratic 
interpolation method, data points for the impedance matrix elements at frequency 
intervals of Afz can be obtained for the entire desired bandwidth. 

The CPU time to evaluate the integral expressions for the voltage vector of a 
uniform current filament, given by Equations (2.57) and (2.58), is similar to the 
CPU time required to evaluate the impedance matrix. Analysis of these voltage 
vector elements with frequency shows that they are slowly varying. For radiation 
problems and loaded antenna computations, this voltage vector is interpolated in 
the same manner as the impedance matrix. 
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2.9.3 Surface Wave Poles 


The integrand for the Zmn integrals are singular when or Dj^ of Equa- 
tions (2.32) and (2.33) equal zero and the integrand for the Vm integrals are singular 
when Dj^ equals zero. The singularities are due to the TE surface waves and 
the Dj^f singularities due to the TM surface waves that exist on the grounded 
dielectric substrate. As the frequency increases, the number of surface wave poles 
can increase. For a lossless dielectric, these poles lie right on the Re(Ac) axis be- 
tween «o, the propagation constant of Region / (in this problem, free space), and 
«2, the propagation constant of Region II. When the dielectric is lossy, the poles 
move into the fourth quadrant, but for low loss dielectrics, they are very close to 
the Re(«) axis. Figure 10 shows a map of the k plane with poles near the Re(/c) 
axis. To avoid these poles without having to compute their location and deforming 
the contour around them, the staggered contour in Figure 10 is used. This method 
has worked well with the imaginary stagger A on the order of 0.1/«o and returning 
to the K axis around 1.1 k2* 

2.9.4 Oscillatory Behavior of the Voltage Vector 

In Subsection 2.5.2, the integral expression for the voltage vector elements 
of a vertical filament of uniform current was derived. The integrand contains an 
exponential function where the argument is a function of (® />?//)• This exponential 
term contributes an oscillatory factor to the integral of Equation (2.59). Figure 11 
is a typical integral with an oscillatory convergence factor. Note that the upper 
limit on this integral is x and the plot is comparing the value of the integral verses 
the upper limit of the integral. This integral oscillates around a converged limit, 
with the size of the oscillations decreasing as x increases. Since the integral of 
Equation (2.59) is evaluated numerically, the infinite upper limit on the k integral 
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must be replaced by a finite limit (which is equivalent to X 2 in Figure 11) where 
the integral is well converged. This limit can be very large, so it would be desirable 
to determine the converged value before the oscillations around this limit become 
negligible. In Figure 11, ®i represents a point where the converged limit of integral 
is apparent. 

The way used to determine this converged limit is by incorporating a peak 
and valley detector in the numeric integration of the Vm integrals. The numeric 
integration computes the integral at discrete points. When a peak or valley is 
found between three integral points, a quadratic interpolation is performed on 
those three points to find the value of the peak or valley. Once a peak and a valley 
are found, the average value of the two is computed. With each new peak or valley, 
a new average is obtained. The integral is considered converged when the percent 
difference between two consecutive averages of a peak and a valley is less than a 
desired tolerance. 
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CHAPTER III 


THE INTEGRAL EQUATION AND MOMENT METHOD 
SOLUTION FOR A VERTICAL THIN WIRE THROUGH A 
GROUNDED DIELECTRIC SUBSTRATE 

3.1 Introduction 

This chapter will describe the integral equation and the moment method [2] 
(MM) solutions to the problem of the self impedance of a vertical thin wire and 
mutual impedance between vertical thin wires on a grounded dielectric substrate. 
The basic geometry of this structure is shown in Figure 12. 

The solution is begun by using the equivalence theorem to replace the perfectly 
conducting thin wire by an unknown surface current J^. Using the thin wire 
approximations, is replaced by a vertical filament of current, [20,21]. Then, 

the integral equation for zl^ is obtained by requiring that the total tangential 
electric field along the filament be equal to zero. Finally, the integral equation is 
solved using the MM. 

3.2 Derivation of the Integral Equation 

This section describes the use of the surface equivalence theorem to derive the 
integral equation for the problem of radiation by a perfectly conducting vertical 
thin wire extending from an infinite perfectly conducting plane and through an 
infinite dielectric layer directly on top of the perfectly conducting plane. The 
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Figure 12: Side view of a vertical, perfectly conducting thin wire extending from 
a perfectly conducting ground plane through a dielectric layer. It is excited by an 
arbitrary source (J\M*) that can be located in either Region / or II. 

geometry of a dielectric layer on top of a perfect conductor will be referred to as a 
dielectric coated ground plane. Figure 12 shows a perfectly conducting thin wire 
of length L in the z direction and radius o (a < A) extending from a grounded 
dielectric substrate. The region above the dielectric layer, referred to as Region /, is 
free space (permeability /iq, permittivity eo). The region inside the dielectric layer, 
referred to as Region //, has a permittivity of 62> permeability of fiQ, and thickness 
T. The origin of the coordinate system is placed at the the point where the vertical 
wire intersects the planar interface of Region / and Region II. The wire is excited 
by source currents (J*,M*) which illuminates the wire with fields (E*,H^). Note 
that (E*, H*) are the fields of (J*,M*) in the presence of the grounded dielectric 
substrate without the wire. (E,H) denote the electric and magnetic fields from 
(J*,M*) in the presence of the grounded dielectric substrate with wire. 




Figure 13: Side view of the equivalent vertical current replacing the perfectly 
conducting thin wire on a dielectric coated ground plane. It is in the presence 
electric and magnetic currents (J^,M*). 

Using the Equivalence theorem and ignoring the currents on the end cap, the 
equivalent surface current on the wire surface will be 

J^ = pxH (3.1) 

Since the wire is thin, surface current on the wire surface is assumed to be z 
directed and uniform with respect to (f>. So J'^ can be written as 


= zl^{z) p' — a 

where the wire current is related to the wire surface current by 


(3.2) 

(3.3) 


The equivalent problem is illustrated in Figure 13. The total electric and magnetic 
fields in either region are 

E = E* + (3.4) 
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H = H' + (3.5) 

where are the fields radiated by the equivalent current in the presence 

of the dielectric coated ground plane. 

On the wire surface, the tangential components of the total electric field must 
equal zero. Then, from Equation (3.4) 

(E‘ + E^) X p = 0 (3.6) 

For the z directed wire this can be written as 

- Ef = E* (3.7) 

Instead of enforcing the boundary condition on the wire surface, it will be enforced 
the center line of the wire. This technique is known as the filamentary test case. 

3.3 The General Moment Method Solution 

This section will use the integral equation developed in the preceding section 
to develop a general MM solution for the geometry of Figure 12. The unknown 
vertical current filament in the integral equation will be replaced by a sum of N 
known expansion functions of unknown strengths. An inner product between each 
side of integral equation and a sum of N known vector weighting functions will 
yield the matrix equation. The expressions for the impedance matrix and voltage 
vector elements will be given. 

To begin the solution for the unknown equivalent vertical current, is rep- 
resented as a sum of N known expansion functions 

= i: (3.8) 

n=l 

where is the n expansion function and In represents the unknown strength of the 
n expansion function. The superscript S is retained to distinguish the expansion 
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function from the unknown strength. Expansion function n is nonzero over the 
range z„i < z < Zn2- Equation (3.7) can be written as 

N 

- E = -E;. (3.9) 

n=l 

is the z directed component of the electric field of expansion function n radi- 
ating in the presence of the grounded dielectric slab. 

By taking the inner product of Equation (3.9) with a set of N known weighting 
functions located on the wire center line, Equation (3.9) is reduce to the matrix 
equation 

(zii/i = m. (3.10) 

where [Z] is an TV x JV impedance matrix, [/] is a TV element current column array 
containing the unknown coefficients /„ of Equation (3.8) and [F] is a iV element 
voltage vector column array. 

The scalar weighting functions employed in this solution will have the form 

= atp = 0 (3.11) 

and thus are chosen identical to the expansion functions. 

Expressions for typical elements of [Z] and [F] are given by 

Z„n = - f’""" Ef„lgdz (3.12) 

and 

Fn»= r^^^EUidz. (3.13) 

Jzml 

3.4 Evaluation of the Impedance Matrix 

In this section, exact integral expressions will be developed for the impedance 
matrix elements using the general expression developed in the preceding section. 


39 



The equations will contain arbitrary expansion functions and arbitrary weight- 
ing functions I^. The electric field of a 4 directed vertical current filament in 
either Region II or Region II is found in plane wave representation. This electric 
field is used in Equation (3.12) to obtain the impedance matrix elements. 

The vector magnetic potential for an infinitesimal vertical current of strength 
I in free space at the origin is 


A = 


Jdz‘ 

^ 

4ir r 


(3.14) 


where 

Ko = wv'eoM (3.15) 

and dz is the infinitesimal length of the current element. The - term in 

Equation (3.14) can be represented as a spectrum of plane waves [20] by 

=L JoWp) (3.16) 

r JQ j'fz 

where Jq is the Bessel function of order 0, the ip in the exponential term is — for 
^ > 0 and -f- for z < 0 and jz Is defined as 


7z = with Re( 72 ) > 0 and Im( 72 ) < 0 (3.17) 


Substituting Equation (3.16) into Equation (3.14) yields 


. Idz' foo 

A = z— 

j47T 


X hm- 


7z 


UdiS 


(3.18) 


If this infinitesimal vertical current element is placed in the dielectric layer of 
the dielectric coated ground plane (Region II) at p' = 0 and z = z\ the vector 
magnetic potential in Region II can be written as, for —T < z < z\ 


A, = — - MM + R_re--"*2(^+T)' 

;4ff h 722 I J 


0dl3 (3.19) 
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and for z' < z < 0, 

Idz' Jq{/3p) 


A,= 


r 

Jo 


jAir Jo 

and for z > 0, 


lz2 


e-nz2{^-^') + R^e^z2^ + il_Te"-''''^2(^+7’) (3,20) 

Idz' f°^ Jo{0p) 


= r-yiLiiToe:>y^^^]f3dfi 

j 4 ir Jo 7z2 


(3.21) 


where, from Equations (3.17) and (3.15), 


7 zl = 

(3.22) 

lz2 = 

(3.23) 

K2 = iVy/e2PQ 

(3.24) 

and, the coefficients Rq,R^'p,Tq are chosen to satisfy the boundary conditions 

E /7 X z = 0 at 2 = —T 

(3.25) 

E /7 xz = E 7 xn at 2 = 0 

(3.26) 

H 77 X z = H 7 X n at 2 = 0 

(3.27) 

where the subscript denotes the region in which the field is 
fields of the vector magnetic potential can be found from 

defined. The E and H 

H = V X A 

(3.28) 

E = — (v2 + k‘^) a 

jive ^ ' 

(3.29) 


where (e, /c) are (eo,Ko) Region I and (e 2 ,« 2 ) for Region II. Substituting 
Equations (3.19) - (3.21) into Equations (3.28) - (3.29) and applying the boundary 
conditions defined by Equations (3.25) - (3.27), the vector magnetic potential of 
a vertical current of infinitesimal length located in Region II becomes, for —T < 
z < z' , 


A, = 


-Idz' foo Jq{/3p) 


r + Xm cos 7^2(2' + T) 

^0 7z2 I J 


j27T Jo 7^2 
e“J7z27’ (.Qg ^^ 2(2 + T) /3 d(3 


(3.30) 
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and for z' < z < {i, 

r +^M<=osj,2{z+T)] 

j27t jo nfz^ '■ ^ 


and for z > 0, 


Az = 


j2ir JO 7^2 
e“JT 22 ^ cos 'yz2{^' + T) fid/3 


-Idz' roo 


(3.31) 


where 


j2rc Jq 

Dm 


r M0P)-^ oosy,2{z‘ + T)e-r<’i‘ffdl3 (3.32) 

JO Dm 


€0722 - €27z1 


Dm 

= C27z1 cos 7^2T + jeo7^2 sm'ifz2T. 


(3.33) 

(3.34) 


If the vertical current of infinitesimal length is located in Region / at a point 
(p = 0 , 2^ = 2^), the vector magnetic potential can be written as, for — T < 2: < 0 , 

^ ^ Jo° + R'_,j,e-^yz2i^+T)j fS dfi (3.35) 


and for 0 < 2 < 2^ 


Az = 


Idz' foo Jq{/3p) 


L 


j4ir Jo 7 zi 


and for z' < z, 


Az = 


Idz' foo Joifip) 


L 




p! -rfziz 


fidfi 


(3.36) 


fidfi 


(3.37) 


j47t Jo 7zi 

where the coefficients R'_>j^,R'q, and Tq are solved for by substituting Equations 
(3.35) - (3.37) into Equations (3.28) - (3.29) and applying the boundary conditions 
of Equations (3.25) - (3.27). The vector magnetic potential for a infinitesimal 
vertical current element in Region I is, for — T < 2 < 0, 


Az = 


Idz' fOO 
j2n 


roo fn t 

L Joifip)YT^-^^^^" cosyz2{z + T) fidfi (3.38) 

JO Dm 
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( 3 . 39 ) 


and for 0 < z < z' , 


Az 


Idz' Jq(^p) 

L n~ 1 cos 7^2^ cos 7^12- 

j2n Jo ifziDM 

eo7z2 sin722r sin7^2«] /3 dfi 


and for z' < z, 


Az = 


Idz' f<x Jq{ 0 p) f ^ 

-r- L p;— C27zl cos7^2^ cos7^i 

j 2 ir Jo IziDm ^ 

^0lz2 sin 722 ^ sin 7^2-2^] d( 3 . 


( 3 . 40 ) 


In general, the vector magnetic potential of an infinitesimal vertical current 
element located in either Region / or Region II can be written as 

Tdz^ /»oo 

= H(^^p)f{(^y^')9{0,z)d/3. ( 3 . 41 ) 

The functions / and g are harmonic functions of z' and 2 respectively. Their 
second derivative with respect to z' or 2 is 


dz'“^ 


f = 


92 

“ 



( 3 . 42 ) 

( 3 . 43 ) 


where 7 = 7^1 for / or in Region I and 7 = 7^2 for / or ^ in Region II. 

Now that the vector magnetic potential for an infinitesimal vertical current 
element is known, the vector magnetic potential for a current of finite length can be 
found. If a filamentary vertical current I{z') exists between z\ and 22, the vector 
magnetic potential is found by integrating the product of I{z') and the appropriate 
vector magnetic potential for an infinitesimal current element with respect to z' 
between z\ and 22. Although the current can be continuous through the Region / 
- Region II interface, it will be analyzed in terms of the sections that wholly exist 
in either Region I or Region II. The expressions generated for the vector magnetic 


43 



potential of the vertical current will have both and 22 W® in either Region / or 
Region II, 

The general form of the vector magnetic potential of a vertical filamentary 
current is given by Equation (3.41). The form of the expression depends on the 
location of 2 . For 2 < 2j < 22 , the general equation can be written as 

A, = ^ h(Mgi{0,z) r Hz’)h(0,z')dz' d0 (3.44) 

and for 21 < 22 < 2 , 

1 /oo /zo 

and for z\ < z < z^, the general equation is written as 

gui(3,z) r I{z')fu{/3,z')dz' 

Jzi 

+ u(0. z) Hz')h(0, z')dz'j d0. (3.46) 

where the L subscript indicates that the / and g for 2 < 2 Ms used and the 
subscript U indicates that the / and g for 2 > 2 Ms used. With the restriction 
that 2 j and 22 are both in the same region and by examining Equations (3.30) - 
(3.31) and (3.39) - (3.40) it is seen that the following relationship exists between 
fu^9Uyfli and gi. 

= 9l{0,z’) = g((3,z') (3.47) 

9U{0>^) = fMz) = fi^,z) . (3.48) 

These relationships are a result of reciprocity and allow Equations (3.44) - (3.46) 
to be written as 

A. = ^ h(0,p)g(0,z) n I(z')mz')dz' d0 (3.49) 


Af — 


1 ro® 
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(3.50) 


A. = j^£°h(/},p)mz) r l(z')g(ff,J)dz' d? 

A. = ^I^h{l3,p)\f(p,z)j\z')g(l3,z')dz' 

+ g{fi,z) I(z’)mz')dz’\dfi. (3.51) 


In order to find the impedance matrix elements for the thin vertical wire, the z 
directed electric field is needed. From Equation (3.29), the Ez field from z directed 
vector magnetic potential is 


Ez 


1 

jwe 



Az. 


(3.52) 


Using Equations (3.49) - (3.52), the general expressions for the Ez fields can 


be written. For z < z\ < Z 2 , 


1 fOO fZo 

= HP,p) 9(0,^)J IU')Hff,z')dz'0Ul3 (3.53) 

and z\ < Z 2 ^ z 


Ez = ^J^°°H^,p)m,z)ri{z')g(0,z')dz'fdl3 (3.54) 

and for z\ < z < Z 2 , 

Ez = ■^J^°°Hl3,p) 0^(^mz)J‘^I(z')g(0,z’)dz' 

+ 9(0 , ») ^ l(^')f(0, ^') <<^') + {^^^g~9(0, 

- ^S^^f(0,z)jl(z)]0d0. (3.55) 

By examining Equations (3.30) - (3.32), the appropriate h,f, and g can be 
determined to write the Ezn held equations for a vertical filamentary wire located 
in Region II with mode current I^ and end points and Zn2~ For —T < z < 

= 2^ /o°° '“7,2(4 + r)F2(4„l,4„2,/|)/j’ d0 (3.56) 
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and for Zni < z < Zn2, 

E,n = r (t-™“ + Xm cos y,2(‘ + T)) 

2 TViV €2 Jo Tz2 ^ 

G2(^nlr«>-^n) + cos 7^2(2 + r)F2(^, «n2) ) 



- }lz2ln] 0 

(3.57) 

and for 

2nl < 2„2 < 2 < 



Ezn - + cos 7 , 2 ( 2 +T)] 

27ro;e2 JO 7z2 ^ 



^ 2 ( 2^111 > 2 ^n 2 j 

(3.58) 

and for 

2 > 0, 



27ro; Jo xJm 

(3.59) 

where 

F2{xi,X2,i) = [e"-^'>'^2® +Xjif cos7z2(® + r)] i(®)d® 

(3.60) 


G2(®1»®2>0 = / cos 7 z 2 (® 4- r)i(®) d®. 

(3.61) 


®1 


From Equations (3.38) - (3.40), the h, /, and g for a vertical filamentery current 
in Region I can be determined and the Ez can be written for a mode current 
with endpoints at z^i and 2 ^ 2 - ^ 2 < 0, 

= r^^^co»7.2(2+r)j'l(s»l,s„2,/»)^*<i/S (3.62) 

27ru? Jo Dm 

and for 0 < z < Znl < - 2 :^ 2 ? 

-1 /•«> Joi^p) f m 

Ezn = 7 , i pr- (e27zlcos722Tcos7^iz 

27TU>eo JO ' 1 z\Dm 

- C07z2sin7227’sin72i2]F’i(2„i,2n25-f^^)/?^<^/3 (3.63) 
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and for Zni < z < z„2> 




-1 r MM (a2^-j 

m 

+ [«2Tzl cos yz2l' cos y^iz - €q7z2 sin jz 2 T sin y^iz] 

h{^,Zn2jS) - nzl^M^ni^)) 1^ 

and for z^i < Zn 2 < z, 

— 1 ’^o(/^P) _-mizA /_ ^ _ _ tS\/3^ 


Ern = 


2iru)tQ Ja 7 ziDm 


where 


Fi(®1,®2>*) = f ^c dx 

Jxi 

- /“®2 

Gi(xi,X2,i) = / [e27zlcos7i2rco8 7*is 


C07z2sin722r sin72ix] t(®) dx. 


(3.64) 


r — (3.65) 

JO '\fAJJM 


(3.66) 


(3.67) 


With the integral expressions for the Em field of a vertical filamentary modes 
in either Region I or Region II, the impedance matrix elements can be obtained 
using Equation (3.12). The self impedance of a thin wire of radius a will be modeled 
as the mutuzd impedance between two filaments of equal length separated by p = a. 
Therefore, the expressions for the impedance matrix elements will be given for two 
modes separated in the p dimension by the distance a. In order to simplify the 
expressions for the elements of the impedance matrix, an I^ mode will be required 
to exist entirely in either Region I or Region II. This is the same requirement 
placed on an I^ mode in the development of the Ezn expressions. The impedance 
matrix element expressions are given for the different locations that an I^ mode 
and an mode can have in the z dimension. For an n mode and an m mode 
both located in Region II with z^n2 ^ ^nl> 

Zmn = /“ — (3-68) 

2it(j}€2 jo yz2 
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and for z„i = and z^2 — ^m2i both the n and m modes in Region //, 


Z 


mn 


= r r*"" [^2 (g-J7z2^ + Xm cos7,2(^ + T)) 

2iru)e2 JO 7z2 J^mi 

<? 2 (*nli*.^n) + / 5 ^' 0 S')'* 2 {* + T)F 2 (z, z„ 2 , li) 

- nz2li\lS,(^)d^l3dl3 (3.69) 


and for z^2 ^ ^ml) the n and m modes in Region II, 

= /'^^^^e-J^*2T^2(^ml,^m2,4)G2(2nl,Zn2,/n)/3^‘^/9 (3-70) 

27 ro;e 2 JQ 7^2 

and for mode n in Region II and mode m in Region /, 

Zm„ = ^ r^^Fl(‘„UZ„2,li)G2(.^nl,^n2,I^)0^d0. (3.71) 

ZTTu; ^0 Vm 


Dm 

For mode n in Region I and mode m in Region //, 


'mn 


= -L/ 

27 ro; •/O 


“ F,(z„i, z„2. )02(*ml, *m2, /^) |9’ <*13 


27ro; •/o 

and for both mode n and mode m in Region I with Zm2 ^ ^nh 


^mn — 


* ■'**'^“’-Fl(j„l,z„ 2 ,/|)< 5 l( 2 ml.*m 2 . 4 )/ 3’<'/3 


27 ro;eo Jo JziDm 


and for z^l = < 2^711 and Zf^2 — ^n2f both modes n and m in Region I, 

Zm„ = 5 -^/”^^^/'"'{/ 3 ^«-"*<^G,(^„,.^,/*) 

27ro;eo JO JzIDm •'^ml 

+ [e 27 zl COS 722 T cos 7 ^ 1 z - eo 7 z 2 sin 722 T sin 72 i 2 ] 


hiz^ZnlJn) - JlzlDMlni^)) ^mi^) dz (3 d(3 


(3.72) 


(3.73) 


(3.74) 


and for z^2 £ 2 ;^!) both modes n and m in Region I, 

Zmn = r ^^^Gl(z„UZn 2 ,lS)Fl{^mU^n. 2 ,Im)ld^ dH. (3.76) 

27ro;eo '0 Tzl^^M 

The impedance matrix elements of Equations (3.68) - (3.75) can also describe 
the mutual impedance between two filamentary vertical currents. If the distance 
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between the two filamentary currents is pQ, replacing the radius a with pQ in the 
preceding impedance matrix element expressions will define the impedance matrix 
element expressions for the mutual impedance. 

3.5 Evaluation for the Voltage Vector 

The wire will be excited by the so-called delta gap generator [19] at z = —T. 
Using this model, the incident electric field of a Fq volt delta gap generator is given 

by 

E* = iVo6(z + T). (3.76) 

Substituting Equation (3.76) into Equation (3.13) results in the expression for the 
voltage vector elements. 

V'm = Vo4(-r) (3.77) 

Equation (3.77) shows that Vm will be nonzero only if the weighting function 
is nonzero at z = — T. 

3.6 Expansion and Weighting Functions 

Now that the general equations of the impedance matrix and voltage vector 
elements are defined, the expansion functions for and the weighting functions 
for will be defined. The wire will be broken into N segments called monopoles. 
Figure 14a shows a vertical wire divided into four segments. There will be N 
modes on the wire, all but one consisting of two monopoles to form a dipole. 
The monopole directly over the ground plane is one mode since its image in the 
ground plane makes it an equivalent dipole. Figure 14b shows the dipole modes 
formed by the four monopole segments on the wire. The functions F and G will 
be evaluated for each region using an arbitrary monopole. The expressions for the 
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electric field at a point between the endpoints of a monopole and the impedance 
matrix elements for overlapping monopoles will be evaluated. 

The vertical wire is broken into a N segments. Figure 14a shows a vertical 
wire on a dielectric coated ground plane broken into four segments. These segments 
have a piecewise sinusoidal function associated with them and are referred to as 
monopoles. The piecewise sinusoidal function is described by 


or 


where 


il{z) 


sin»Crf(z2-2) 

sin 


for zi < z < Z 2 


0 elsewhere 


*2(2) 


sin ) 

sin/c^Az 

0 


for z\ < z < Z 2 
elsewhere 


= 

Az = 
k^Az < 


KQ for Region / 
K 2 for Region II 

Z2 - zi 
IT 
2 ' 


(3.78) 

(3.79) 

(3.80) 

(3.81) 

(3.82) 


The F and G functions shown in the integrands of Equations (3.56), (3.58), 
(3.59), (3.62), (3.63), and (3.65) are obtained for the monopoles of Equations (3.78) 
and (3.79) between points z\ and Z 2 - For the function, 

^ — ■JT' - [-J'izl sin kqAz + kq cos kqAz] (3.83) 

P^sxiikqAz ' ' 

^ (c27z 2 cos 7i2r[/co (cos 7^122 - cos 7^121 cos /cqAz) 
sin KQ Az 

+ 7ai sin 7^121 sin «o^2] 

+ eo7z2sin7i2T[Ko(sin72iZ2 - sin72izi cos «oAz) 
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_ ^ ^ ^ ^ ^ y ^ y 


a) A vertical wire through a grounded dielectric substrate broken into four seg- 
ments. 



b) The piecewise sinusoidal functions on each segment to construct the MM modes. 

Figure 14: An example of a MM modal model for a vertical wire through a 

grounded dielectric substrate. 
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( 3 . 84 ) 


- 72 l cos 72 l^t sinfcoAzj) 

h = -:v 5 ^ («2 [e"-''>'^222 - cos/C2A2e“-'T'*2^l| 

^ /32 sin k 2 Az V ^ I ^ J 

— jjz 2 sin K 2 Az€~^^^^^ 

+ [«2( cos 7.22 (22 + T) - cos cos 7^2(21 + T)) 

+ 722sin/C2Azsin7^2(-^l + r)l) (3.85) 

^2 = . — — ^ (k 2 [cos 7^2(22 + ^) -cos K2 Az cos 7^2(21 +T)] 

)0'^sin/C2Az 

+ 7^2 sin K2A2 sin 722(21 + T)) ( 3 . 86 ) 


and for the »2 segment function, 

-^1 == a? • ^ — A~ - [j7zl sin/CflAz + /C() cosacqAz] ( 3 . 87 ) 

P'^sin/CQAz ' ' 

= :b2 2 . " — r-(e 27^1 cos 722 T’[«o(cos 72 izi - cos 72 iZ 2 CosKoAz) 

Sin kq Az 

- 72I sin 72 iZ 2 sinKo^ 2 ] 

+ co7z2sin7227’[«o(sin72izi - sin72iZ2 cos / cqAz) 


F 2 = 


+ 721 cos 72 iZ 2 sinK 0 ^ 2 ]) 

1 


( 3 . 88 ) 


^«2 |e — cos«2Aze ^"^z2^2 


13 ^ sin /C2 Az 
+ J 722 sin «2 Aze"'^'*'22^2 

+ -^Af ['*2 (cos 722(21 + T) - cos /C2 A z cos 722 (22 + T)) 

- 722 sin «2 Az sin 722(22 + T)]) 

1 


( 3 . 89 ) 


G2 = («2 [cos 722(21 +T)- cos K2 Az cos 722(22 + r)] 

- 722 sin K2 Az sin 722(22 + T)) . ( 3 . 90 ) 


When a dipole is constructed from two monopoles, some terms in the F and G 
will cancel. These terms are due to the discontinuous current in the monopole. The 
first step for identifying terms that cancel is to determine the z directed electric 
field of a dipole. For this step the terms that cancel are called charge terms, since 
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the physical representation is a charge at the end of the monopole where the current 
is non zero. The integrand of the spectral integral for the Ez field of a z directed 
monopole consists of terms that have & I3~^ dependance and a dependance. 
When two monopoles are combined to form a dipole, all 0^ dependent terms are 
charge terms and cancel regardless of the position of the monopoles or the field 
point. Without removing the charge terms, the integrals of the monopole are non 
convergent. All of the charge terms are products with either 7^1 and sin/coAz in 
them for Region I or 7^2 sin«2Az in them for Region II. As an example, for 
Fi of the ii function, the charge term is —j^zi sinKoA2c”-^*i^i when Fi is in the 
integral equation for the i directed electric field. 

The last step in identifying terms that cancel requires the determination of 
the dipole to dipole impedance matrix element integral expressions. The terms 
that cancel from these expressions are of similar algebraic form as charge terms, 
and are called current discontinuity terms. Although the integrals will converge 
with the current discontinuity terms, removing them can decrease the convergence 
time of the numerical integration. The identification of the current discontinu- 
ity terms involves examining all possible dipole to dipole configurations. These 
configurations are dependent upon the z location of the monopoles. 

In Equations (3.57) and (3.64), the F and G are functions of the z component 
of the field point. Evaluating F and G for the piecewise sinusoidal expansion 
function and expanding the part of the integrands in Equations (3.57) and (3.64) 
which are a function of z, the integrands can be simplified. For (3.57), the part 
of the integrand that contains the functions of z is 


02 (e-rrz2^+XMCosjz2i^ + T)) G2(z„i,z,/|) 
+0^ cos 722(2 + T)F 2 {z,Zn 2 ,In) - 31z2ln 
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(3.94) 


Letting = i\, evaluating F and G, and simplifying the result, Equation (3.91) 
becomes 

siniz^z [7^2sin7^2(^n2 + ^) ~ «2 cos 7^2(2,11 + r)cosK2^z] 

005722(2 + T) K2e“^'>'«22n2 -I- Xjif (k 2 COS722(2n2 + r) • (3-92) 

+ 7^2 sin Jz2{^nl + T) sin K 2 AZ - K 2 cos Jz2(^nl + T) cos K2^-2)]) 

Doing the same with = * 2 > Equation (3.91) becomes 

slnigAi (^"•''’'^22 [«2 cos7^2(2nl + T)] + 0037^2(2 + T) [(-J722 sin «2 Az 
-«2 cosk2 ^2) + XM{-'yz2 ^^^Iz2{^n2 + T) sin/C 2 Az 

-«2 cos7^2(2n2 + T)cosk 2 Az + K2 cos7^2(2„i + T))]) 

(3.93) 

For Equation (3.64), the part of the integrand which contains functions of z is 
/?^e-JT'^i^Gi(z„i,z,/|) + /?2 [627^1 cos 7^2 ^’cos 7^1Z 
-e07z2 smjz 2 T sin 7^1 z] Fi(z, z„ 2 ,In ) “ J7zl^A/-^n 
Letting = i\ and evaluating, Equation (3.94) becomes 

sin/cpAz [C27zl cos'fz2T {-kq cos-f^iZni cos/cqAz 

+721 sin72iz„i sinKoAz) - £9722 sin 722^ {-kq sin72iz„i cos /cqAz 
- 721 cos72iz„i sin kqAz)] + [e27^l COS72IZ 
-60722 sin72iz] 

Letting = 12 , Equation (3.94) becomes 

slnipAz [e27zlcos722Tcos72iz„i - 69722 sin722rsin72iz„i] 

- [62721 cos 722 ^ cos 72 IZ - 60722 sin 722r sin 72iz] /coe“-^'>' 2 i-^»i 2 cos ko Az • 
sin /coAzj j 

(3.96) 

Note that the charge terms in the above expressions are the terms that contain 
72 I and sin^oAz or or 722 and sin/C 2 Az. 


(3.95) 
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By using Equations (3.92) - (3.93) in Equation (3.57) and applying this result 
to Equation (3.12), the integral expression for an impedance matrix element of 
overlapping piecewise sinusoidal monopoles with expansion mode and weighting 
mode can be written. For the monopoles in Region II and with I^=i\, 

" 2^1 7 z 2 (Fo{zmUZm2, 4) bz2 sm7z2(^nl + T) 

- K2COS722(«nl + r)cosK2Az] + G2{zml,Zm2Jm) [/C2e“-^*2^n2 

+ («2COS722(«n2 + ^) + 7z2 sin7z2(«nl + T)s\nK2Az 

- K2 cos ‘yz2{znl + T) cos K2 Az)]) fid/S. (3.97) 

For the monopoles in Region II with /n=i2» 

Zmn = (^0(2ml, 2,„2,7 ^) [k 2 COS7z2(«n2 + T’)] 

27ra;c2 *'0 7z2 sm K2 Az '' 

+ G2{zml,Zm2Jm) [(-J7z2 sin K2Az - K2 cosK2Az) 

+ («2 cos7z2(znl + T) - 7^2 sin7z2(zn2 + T) sin K 2 Az 

- K2 C0S7z2(«n2 + T) cos /C2 Az)]) d/3 (3.98) 

where 

Fo(xi,« 2 j*)= f ^ e~^^^^i(x)dx. (3.99) 

Jxi 

Note that is the same as Fi with 7^1 replaced by 7z2* For a piecewise si- 
nusoidal monopole, the integral in Equation (3.99) can be obtained from Equa- 
tions (3.83) and (3.87) by replacing 7^1 by 7^2 and kq by /C2* Using Equations (3.64) 
and (3.95) - (3.96) the impedance matrix element for overlapping monopoles in Re- 
gion / becomes, with 

„ 1 /“oo Jn(l3a) ,, ^/ . .A 

Zmn = 7i L : ( eo7z2 cos7z2J (7^1 sm7ziZni sin kqAz 

27tu;co jo 7zi sm kq Az 

- KQ cos 7 ziz„i cos kqAz) -I- co7z2 sin7z2r («o sin 7 ziz„i cos kqAz 
+ 7zl C 0 S 7 ziz„i sin/coAz)] Fi{zml,Zm2Jm) 
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(3.100) 


For the monopoles in Region / with /^=*2> 

Zmn = TT- — L (A(^ml» ^m2> [^07z2 cos722^ cosTzHnl 

27ra;eo ./O 72ismKo^'®^ 

- eo7z2sin7^2r sin 72 iz„i] -f Gi(zrnl> 2m2> ^m) cos kqAz 

+ sinKgAzj^ /3rf/3. (3.101) 

3.7 Numerical Difficulties 

The integrand for the Zmn integrals are singular when Dj^ equals zero. The 
Djlf singularities is due to the TM surface waves that exist on the dielectric coated 
ground plane. As the frequency increases, the number of surface wave poles can 
increase. For a lossless dielectric, these poles lie right on the Re(>c) axis between 
«0) the propagation constant of Region I (in this problem, free space), and « 2 , the 
propagation constant of Region II. When the dielectric is lossy, the poles move 
into the fourth quadrant, but for low loss dielectrics, they are very close to the 
Re(«) axis. Figure 10 in Subsection 2.9.3 shows a map of the k plane near the 
Re(«) axis, including a staggered contour which avoids the surface wave poles. 
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CHAPTER IV 


NUMERICAL RESULTS 


4.1 Introduction 

This chapter compares calculated and measured input impedance and radar 
cross section (RCS) of unloaded and loaded rectangular patch microstrip anten- 
nas. The mutual impedance of vertical wires in a grounded dielectric substrate is 
compared to previously calculated data [21]. The first section will compare calcu- 
lations to a compact range measurement of the RCS of a rectangular patch on a 
grounded dielectric substrate. A calculation will be made for input impedance of 
this patch when corner fed in order to show the modal relationship between RCS 
peaks and input impedance resonances. The patch will then be loaded near the 
corner with a 50 ohm load. The RCS measurement and calculation for the loaded 
microstrip patch will be compared. The next section will compare the calcula- 
tion to the measurement for the input impedance of a microstrip antenna. This 
will include a wide band comparison and near resonance narrow band compar- 
isons. The antenna will be loaded with a 50 ohm load near the corner and the 
input impedance calculation and measurement compared. Next, the calculation 
of the mutual impedance between two vertical thin wires imbedded in a grounded 
dielectric substrate will be compared to previously calculated data [21]. 
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Figure 15: The rectangular microstrip antenna used for RCS calculation and 

measurement. 

4.2 Plane Wave Scattering Results 

This section will compare the MM computation for the radar cross section 
(RCS) of a microstrip patch to measurement. Figure 15 shows the geometry of the 
microstrip patch used for the unloaded RCS data. Figure 16 shows the geometry 
of the same patch with location of the vertical filament of current used for an input 
impedance calculation. Figure 17 shows the same patch used for loaded RCS data. 
The patch is loaded with 50 ohms for a comparison of calculated and measured 
RCS data. The measured patch was loaded through a coeucial feed. 

The patch has dimensions of 2H = 3.66 cm and 2W — 2.60 cm. The substrate 
has thickness T = 0.158 cm, relative permittivity Cr = 2.17, and a loss tangent of 
0.001. The incident plane wave has an angle of incidence of — (60°, 45®) 

and its electric field is 6 polarized. Figure 18 shows the RCS of the patch from 
2 to 10 GHz. Figures 19 and 20 show the calculated real and imaginary parts of 
the input impedance of this patch over the same frequency range when excited by 
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f Yf ) = (1 • 67 cm, 1.14 cm) 


Figure 16: The rectangular microstrip antenna of Figure 15 with a feed port for 

input impedance calculation. 
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={ 1^33 cm, 0.80 cm) 


Figure 17; The rectangular microstrip antenna of Figure 15 loaded for RCS 

calculation and measurement. 
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a vertical current located at (®y,y^)=(1.60 cm, 1.10 cm), which is very close to a 
corner of the patch. Figure 21 compares the wideband RCS of the same patch 
when loaded by 50 ohms at (®i,2/i)=(1.33 cm, 0.80 cm). 

The MM computations of Figures 18-21 were made with A/c = 400 MHz 
(the frequency interval for impedance matrix and wire voltage vector computation) 
and Afz = 5 MHz (the frequency interval for the MM solution) [14]. The solid 
curves are the calculated data with the dots on Figure 18 representing the points 
where the impedance matrix has been calculated. The number of modes used to 
calculate the data of Figures 18-20 were 17 for the frequency range of 2 to 6 GHz 
and 31 from 6 to 10 GHz. Note that the changing of the number of modes causes 
a discontinuity at 6 GHZ for some of the calculated data. To test for convergence, 
the MM solution was run for 84 modes. The data was virtually the same. The 
total CPU time to calculate the data in Figure 18 was about 27 minutes on a VAX 
11/780. The computation of the data in Figures 19 - 20 took about 30 minutes 
on a VAX 8550. For Figure 21, 4 modes were used to calculate the curve from 2 
to 5 GHz, 12 modes from 5 to 7.5 GHz, and 24 modes from 7.5 GHz to 10 GHz. 
Discontinuities in the calculated data at 5 GHz and 7.5 GHz are due to changing 
the number of modes. The total CPU time to compute the data of Figure 21 was 
about 30 minutes on a VAX 8550. 

Comparing Figures 18 - 20 shows that the same modes can be excited by both 
sources. The edge feed location was chosen to excite all modes in the displayed 
frequency range. Some of the peaks in the RCS and resonances in the input 
impedance are caused by a resonance of a single patch mode, others are caused by 
the mutual coupling between patch modes. This will be analyzed in more detail 
in Section 4.3. 
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Figure 18 : A comparison of the computed and measured theta polarized RCS of 
a rectangular patch on a grounded dielectric substrate. 
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2H=3.66cm 2W=2 . 60cm 
£= 2.11 tan(5=0.001 T=0. 158cm 


Xj=1.67cm y^=1.14cm 



F (GHz) 


Figure 19: A calculation of the real part of the input impedance for the 
microstrip antenna of Figure 16. 
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2H=3 . 66cm 2W=2 . 6cm 


€^=2.17 tan6=0.001 T=0. 158cm 


Xf=1.67cm y^=1.14cm 



Figure 20: A calculation of the imaginary part of the input impedance for the 

microstrip antenna of Figure 16. 
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2H = 3.66cm 2W = 2.60cm 
Gj. = 2.17 tan(5 = 0.001 T = 0.158cm 
= 1.33cm = 0.80cm = 50 Q 



F (GHz) 


Figure 21: A comparison of the computed and measured theta polarized RCS of 
a loaded rectangular patch on a grounded dielectric substrate. 
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4.3 Input Impedance Results 


This section will compare the calculated input impedance of a microstrip an- 
tenna excited by a filament of uniform current to the input impedance of a coaxially 
fed microstrip antenna. The measured antenna is loaded through a coaxial feed. 
Note that the current filament is not a model for the coaxial feed since it does 
not account for the aperture in the ground plane and the current on the outer 
conductor of the coax. Figure 22 shows the geometry of the unloaded microstrip 
antenna used to generate Figures 25 to 31. Figure 23 shows the geometry of the 
edge loaded microstrip antenna used to generate Figures 32 and 33. Figure 24 
shows the geometry of the corner loaded microstrip antenna used to generate Fig- 
ures 34 and 35. The thickness of the substrate was measured to be 0.088 cm. 
The loaded antenna had a load impedance of 50 ohms. The relative permittivity 
and loss tangent, as supplied with the substrate, are 2.17 and 0.0015 respectively. 
Referring to Figure 5a, the feed is located at = (0.85 cm,1.22 cm). Fig- 

ure 25 shows the wide band input resistance and Figure 26 shows the wide baud 
input reactance of this microstrip antenna. Figures 27-31 show the calculated 
and measured narrow band input impedance plotted on a Smith Chart for each 
resonance between 1 and 7 GHz. Figure 34 shows loaded input impedance from 
1 to 4 GHz with the load located at {xi,yi)={1.53 cm,l.l cm). Figure 35 shows 
the loaded input reactance. 

For the calculated data of Figures 25 and 26, A/c was chosen to be 400 MHz 
and Afz was chosen to be 5 MHz. Twelve modes were used in the MM solution to 
compute the curve between 1 and 3.5 GHz, 24 modes were used for the calculation 
between 3.5 and 6.4 GHz, and 50 modes were used for the calculation of the curve 
between 6.4 and 8.0 GHz. The convergence of this data was tested by using 24 
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= (0.85 cm, 1.22 cm) 


Figure 22: The rectangular microstrip antenna used for input impedance 

calculation and measurement. 
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(^ftYf) = (0.85 cm, 1.22 cm) 

(x^rY^) = (0.0 cm, -2. 4 cm) 

Figure 23: The rectangular microstrip antenna of Figure 22 loaded at the edge 
for input impedance calculation and measurement. 
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(Xf/Yf) = (0.85 cm, 1.22 cm) 

(^l/Yx.) “ (“1.5 cm, “2.2 cm) 

Figure 24: The rectangular microstrip antenna of Figure 22 loaded at the corner 
for input impedance calculation and measurement. 

modes in the range from 1 to 3.5 GHz, using 32 modes from 3.5 to 6.4 GHz and 
using 72 modes from 6.4 to 8.0 GHz. The curves were virtually the same after the 
increase in modes. The total CPU time needed to generate the calculated data of 
Figures 25 and 26 was about 37 minutes on a VAX 8550. 

For the calculated narrow band data of Figures 27 - 31, A/c was chosen to 
be 200 MHz and Afz was chosen to be 1 MHz. The number of modes used to 
compute the data of Figures 27 and 28 was 4. For Figure 29, 12 modes were used, 
for Figure 30, 32 modes were used and for Figure 31, 72 modes were used. The dots 
on the calculated and measured curves represent the measured frequency closest 
to resonance, and the frequencies 1% and 2% above and below the closest resonant 
frequency. 

The calculation for loaded input impedance of Figures 34 and 35 was made 


(Xf ,Yf )• 


5.0 < 

cm 

♦-3.4 cm 


•(x^rYx.) 
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2H=3.4cm 2W=5.0cm 
€^=2.17 tan(5=0.0015 T=0. 08779cm 
Xj=0.85cm y^=1.22cm 



F (GHz) 


Figure 25: A comparison of the computed and measured real part of the input 

impedance for a microstrip antenna. 
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2H=3 . 4cm 2W=5 . 0cm 

£^.=2.17 tan<5=0.0015 T=0. 08779cm 

x_=0.85cm y =1.22cm 
f ^ f 



Figure 26: A comparison of the computed and measured imaginary part of the 
input impedance for a microstrip antenna. 


2H=3.4cm 2W=5.0cm 
er=2.17 tan(5“0.0015 T=0. 08779cm 
Xj=0.85cm y^=1.22cm 



- MEASURED 

- CALCULATED 


Figure 27: A comparison of the computed and measured input impedance around 
first resonance plotted on a Smith Chart. 
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2H=3.4cm 2W=5,0cm 
er=2.17 tan6=0.0015 T=0. 08779cm 
x^=0.85cm y^*1.22cm 



MEASURED 

CALCULATED 


Figure 28: A comparison of the computed and measured input impedance around 
second resonance plotted on a Smith Chart. 
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2H=3.4cm 2W=5.0cm 
e^=2.17 tan6=0.0015 T=0. 08779cm 
x^=0.85cm yj=1.22cm 



- MEASURED 

- CALCULATED 


Figure 29: A comparison of the computed and measured input impedance around 
third resonance plotted on a Smith Chart. 
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2H=3.4cm 2W=5.0 

€r=2.17 tan6=0.0015 T=0. 08779cm 
x^=0.85cm y^=1.22cm 



MEASURED 


- CALCULATED 


Figure 30: A comparison of the computed and measured input impedance around 
fourth resonance plotted on a Smith Chart. 
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2H=3.4cm 2W=5.0cm 


Gr=2.17 tan6=0.0015 T=0. 08779cm 


x^=0.85cm y =1.22cm x^=0.0cm y^=2.4cm z =50 fl 

f f L L li 



Figure 32: A comparison of the computed and measured real part of the input 
impedance for an edge loaded microstrip antenna. 
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2H=3.4cm 2W=5.0cm 


€^=2.17 tan6=0.0015 T=0. 08779cm 
x-=0.85cm y =l,22cm x^=1.5cm y =2. 2cm z =50 

f ^ f L L L 



Figure 34: A comparison of the computed and measured real part of the input 
impedance for a corner loaded microstrip antenna. 
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2H»3.4cm 2W=5.0cin 


€j.=2.17 tan6=0.0015 T=0. 08779cm 

Xj=0.85cm y^=1.22cm yj_=2.2cm Zj^=50 ^ 





2H = 3.4cm 2W = 5.0cm 

6r = 2.17 tand = 0.0015 T = 0.08779cm 
x^= 0.85cm y^= 1.22cm -1.5cm Yi,- -2.2cm Zj^- OQ 



F (GHz) 


Figure 36; A comparison of the measured real part of the input impedance for an 
unloaded and corner shorted microstrip antenna. 
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2H = 3 . 4cm 2W = 5 . 0cm 

= 2.17 tan6 = 0.0015 T = 0.08779cm 
x^= 0.85cm y^= 1.22cm x^= -1.5cm -2.2cm Zj_= 0(1 



Figure 37: A comparison of the measured imaginary part of the input impedance 
for an unloaded and corner shorted microstrip antenna. 
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2H = 3 . 4cm 2W = 5 . 0cm 

€r = 2.17 tan(5 = 0.0015 T = 0.008779cm 
x^= 0.85cm yf= 1.22cm Xl= -1.5cm yj^= -2.2cm Zj^= Of) 



F (GHz) 


Figure 38: A comparison of the computed and measured real part of the input 
impedance for a corner shorted microstrip antenna. 
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2H = 3.4cm 2W = 5.0cm 


£r = 2.17 tan6 = 0.0015 T = 0.08779cm 
Xf= 0.85cm Yf= 1.22cm x^= -1.5cm yj^== -2.2cm Z^= Ofi 



Figure 39: A comparison of the computed and measured imaginary part of the 
input impedance for a corner shorted microstrip antenna. 
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with A/c = 400 MHz and A/* = 6 MHz. The total number of modes used was 4. 
The CPU time require to generate this data on a VAX 8550 was about 4 minutes. 

Figures 36 and 37 show a comparison of input impedance for the measured 
unloaded and shorted microstrip antenna. The geometry of this antenna is shown 
in Figure 24 where the load port is now a vertical wire from the patch to the 
ground plane. One might expect a short near the corner to eliminate resonances, 
particularly the first three since they were greatly reduced by the impedance load. 
Figures 36 and 37 clearly show that this is not the case. Figures 38 and 39 show 
a comparison of calculated input impedance and measured input impedance for 
the shorted microstrip antenna. The calculation was made using 4 modes from 
1.5 GHz to 2.5 GHz, 12 modes from 2.5 GHz to 4.5 GHz, 17 modes from 4.5 GHz 
to 6.2 GHz, and 32 modes from 6.2 GHz to 7.5 GHz. 

The calculated and measured data shown are close in resonant frequency but 
usually differ in magnitude. This can be accounted for by the following. First, the 
excitation for the calculation is not a model for a coaxial feed. Second, the vertical 
filament of current excitation does not include an attachment mode between the 
filament and the patch. Third, the microstrip antenna has very high Q resonances 
and any loss not modeled in the MM solution could result in lower magnitudes for 
the measurement. When looking at the calculated results, particularly the input 
reactance data, the discontinuities in the curves are the result of adding more 
modes to the MM solution. The higher order modes account for the inductive 
shift of the data as the frequency goes higher. 

It is interesting to examine the relationship between the impedance matrix and 
the input impedance. One might expect an input impedance resonance to directly 
correlate with a resonance of a particular mode or modes in the impedance matrix. 
This is the case for three of the five resonances shown in Figures 27-31. The other 
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two are caused by the mutual coupling between modes. 

For Figure 27, the most strongly excited mode is the 1 0 directed mode and it 
is resonant. The most strongly excited and resonant mode for Figure 28 is the 1 0 x 
directed mode. The resonance shown in Figure 30 has the 3 0 y directed mode most 
strongly excited and resonant. This is no surprise since it is two octaves above the 
resonance caused by the 1 0 y directed mode. The resonance shown in Figure 29 
does not correspond to a resonant mode. The most strongly excited mode is the 
1 1 X directed mode. By analyzing the input impedance where two modes were 
used for calculating the impedance matrix, it was found that the 1 1 ic directed and 
1 1 y directed modes were needed to create an input impedance resonance at that 
frequency. Neither of these modes were near resonance. Similarly, for Figure 31, 
there is no resonant mode in the impedance matrix. The most excited mode is the 
3 1 y directed mode. The modes that are required to show this resonance in the 
MM computation were not determined since a large number of modes were needed 
to generate it for Figures 25, 26, and 31. This resonance was most likely caused by 
the interaction of the 3 1^ directed mode and several other higher order modes. 

4.4 Mutual Impedance Between Vertical Thin Wires in a Grounded 
Dielectric Substrate 

The mutual impedance calculation with respect to their separation for two 
vertical thin wires of equal length inside the dielectric region is compared to previ- 
ously calculated data [21]. The separation and dimensions have all been normalized 
to the free space wavelength at the operating frequency. The thickness of the sub- 
strate is 0.4152 Aq. The dielectric constant of the substrate is 2.45 and the substrate 
is considered lossless. The wires extend from the ground plane at 2 = — 0.4152Aq 
to 2 = — 0.2636Ao and have a radius of O.OOIAq. The separation of the wires is 
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varied from 0.03Ao to 1.95Ao. For the published curve, the excitation was a mag- 
netic frill current around the aperture of the coax. Figures 41 and 42 show that 
the so called delta gap generator yields virtually the same curves for the real and 
imaginary parts of the mutual impedance. The geometry of this calculation is 
shown in Figure 40. 

For the MM computation, the wires were broken into three equal length seg- 
ments. This corresponds to three dipole modes on the wires. The data was com- 
puted at intervals of 0.01 Ao separation. The numeric integration had a maximum 
limit on the k axis of 200kq but the actual maximum value of k occurred when 
the integration converges to within .2%. As the wires moved further apart, the 
maximum value of k for convergence decreases but the density of the integration 
points must increase since the Jq Bessel function in the impedance matrix integrals 
oscillate more rapidly with respect to k as the separation s increases. The total 
CPU time to compute the data in Figures 41 and 42 was about 9 minutes on a 
VAX 8550. 

The mutual impedance is defined from the two port problem. The excitation 
ports are at the base of each wire at z = — 0.4152Ao. The two port linear equation 
is written as 


^11 

Z12 


h 



■^21 

Z22 


h 


V2 


where Z\\ and Z 22 are the self impedances of wire 1 and wire 2 respectively, Z \2 
and Z 21 are the desired mutual impedance, I\^V\ are the current, voltage at port 
1, respectively, and J 2 , V 2 are the current, voltage at port 2, respectively. The two 
port impedance matrix is developed by solving the MM matrix equations first with 
the excitation at port 1 and then with the excitation at port 2. When the MM 
matrix equations are solved with the excitation at port 1, the currents at port 1 
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L=0.1516>\o 

T=0.4152>\o 


Figure 40: Two thin wires imbedded in the dielectric of a grounded dielectric 
substrate. The wires have three dipole MM modes on them. 
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and port 2 can be determined. Dividing these currents by the voltage exciting port 
1, the Yii and Y21 elements of the two port admittance matrix are determined. 
Similarly, the excitation at port 2 will define the Y12 and Y22 elements of the 
two port admittance matrix. By inverting this admittance matrix, the two port 
impedance matrix is obtained. 

Figure 43 was constructed to illustrate the effect of the dielectric substrate 
on the mutual impedance. The dielectric substrate was replaced by free space 
for the same geometry and the mutual impedance was computed over a wider 
range of separation. The solid curve is the mutual impedance for the two wires 
imbedded in the grounded dielectric substrate and the dashed curve is the free 
space mutual impedance. The mutual impedance magnitude falls off as roughly ^ 
when the the wires are in free space. This is as expected since the magnitude of the 
mutual impedance is linearly proportional to the electric field and the electric field 
decreases as ^ in free space. When the dielectric is present, there is substantially 
more coupling between the two wires as the separation increases. The is due to 
the TM surface waves on the dielectric. 
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Figure 42: A comparison of present calculation to published data [21] for the 
imaginary part of the mutual impedance between two thin wires imbedded in a 

grounded dielectric substrate. 
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Figure 43: A comparison of the magnitude of the mutual impedance of Figure 40 
and the mutual impedance of the same geometry where the dielectric substrate is 

replaced by free space. 




CHAPTER V 


SUMMARY 


This report has described an integral equation and MM solutions to the prob- 
lems of scattering and radiation from an unloaded and loaded rectangular micro- 
strip patch and the self and mutual impedance of vertical thin wires protruding 
through a grounded dielectric substrate. 

For the rectangular microstrip patch, an integral equation was developed for 
unknown surface currents on the patch. The MM was applied to represent this 
integral equation as a system of linear equations. A crucial part of the solution 
was the integral expressions for the impedance matrix elements. The integrands 
in these expressions were analyzed for even-odd symmetries in order to reduce 
the time of numeric integration. The voltage vector integral expressions for a 
uniform vertical current were also analyzed for even-odd symmetries to reduce the 
time required for their computation. The far field of each mode on the patch was 
expanded asymptotically so that the voltage vector for an incident plane wave 
and the scattered or radiated far field could be expressed analytically. Coaxially 
loading of the patch was modeled by a lumped load in a vertical uniform current 
filament. The current through the load is determined by a two port analysis. 

This report used another technique to improve computational efficiency of the 
MM solution. The impedance matrix and voltage vector elements were analyzed 
with respect to their frequency dependance. It was found that they varied slowly 
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and smoothly with frequency so an interpolation scheme was employed to generate 
impedance matrix and voltage vector data at frequency intervals in between those 
that were calculated [16]. The result was the ability to generate wide band data 
from relatively few impedance matrix computations. Calculated and measured 
data were compared for wide band RCS, wide band input impedance and narrow 
band input impedance of both unloaded and loaded microstrip antennas. The 
RCS curve showed good agreement between calculation and measurement. It also 
showed how the impedance matrix interpolation could discern the narrow band 
peaks in the RCS that were not obvious from the points where the MM calculated 
the impedance matrix. The data presented for the wide band input impedance 
also showed good agreement between calculations and measurement. 

For the vertical thin wire in the grounded dielectric substrate, an integral 
equation was developed for surface currents on the wire. Since a thin wire has a 
diameter much smaller than a wavelength, this surface current was approximated 
by a z directed filament of current. The current was obtained by the MM. The 
excitation chosen to determine self and mutual impedance was the so called delta 
gap generator. The mutual impedance between two thin wires imbedded in the 
dielectric is compared to a previous calculation [21] which used a magnetic frill 
current for excitation. The close agreement between the calculated and published 
curves show that the two types of excitation yield virtually identical results. 
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